1 


l 





WtAW COPY* RETURN ] 

APWL (WLOE) 

XIRtLAND APB, R MEX 


THE EFFECTS OF A NONRIGIDLY 
SUPPORTED BALLAST ON THE 
DYNAMICS OF A SLENDER BODY 
DESCENDING THROUGH THE ATMOSPHERE 

V> ; » A . ' , 

by Gerald N. Malcolm ; 

Ames Research Center *' 

Moffett Field, Calif. 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • JANUARY 1970 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



□132M55 


1 . 

4. 

7. 

9. 

12 . 


Report No. 

NASA TN D-5622 


2. Government Accession No. 


Title and Subti tie 

THE EFFECTS OF A NONRIGIDLY SUPPORTED BALLAST ON THE DYNAMICS 
OF A SLENDER BODY DESCENDING THROUGH THE ATMOSPHERE 


3. 

5. 

6 . 


Author(s) 8. 

Gerald N- Malcolm 


Performing Organization Name and Address 10- 

NASA Ames Research Center 
Moffett Field, Calif. , 94035 

11 . 


13 . 


Sponsoring Agency Name and Address 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Washington, D. C. , 20546 

14. 


Recipient's Catalog No. 


Report Date 

January 1970 

Performing Organization Code 


Performing Organization Report No. 
A- 3 291 


Work Unit No. 

124-07-02-14-00-21 

Contract or Grant No. 

Type of Report and Period Covered 
TECHNICAL NOTE 


Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

The effects of a nonrigidly supported ballast on the dynamic behavior of a slender body with linear aerodynamics 
descending through the atmosphere in planar motion are investigated. The equations of motion of both the main vehicle and the 
spring- mounted ballast are derived and specific examples are presented to illustrate the resulting motions. Divergent 
oscillations of the main vehicle occurred for those conditions near the coalescence of frequency of the main body and 
ballast motions. Divergence occurred earlier in the flight when spring stiffness was decreased, or when ballast damping 
was increased for a given spring stiffness. 


17 . 


Key Words Suggested by Author 

Dynamics of slender entry body with non rigid ballast 


18. Distribution Statement 

Unclassified - Unlimited 


Security Classif. (of this report) 


20. Security Classif. (of this page) 


21* No. of Pages 


22 . 


Price* 


Unclassified 


Unclassified 


20 


S 


3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 





SYMBOLS 


A 

b 

Cm 

c m 0 

Cm c 

CMQ 


sys 


c m q + C m . 


C N n 


d 

fr 


h 

I 


L sys 


reference area, , ft 2 

damping coefficient for ballast system, slugs/sec 

, . ^ ^ pitching moment 

pitching-moment coefficient, r ^ —g- 


3C, 


pitching-moment curve slope of cone alone (without ballast) , 

pitching-moment curve slope of system (cone plus ballast) 

C N, 


’in 


damping-in-pitch parameter C m + C m . - 


"a 


3C 


m 


m^d 2 / I S y S 

9 CL 


J m 


damping-in-pitch derivative, 3 - ^ d/v y + 3 ( a d/ v) 

normal-force curve slope 
axial-force coefficient 

reference length (diameter of cone base), ft 

natural frequency of ballast attached to platform of infinite 

c l [Y 

mass, f n = — / — , cps 
°y 2tt / m2 r 

natural frequency of cone with rigidly mounted ballast, 

f = _L / q °°^ d a sys ~ 

in o / -r > L 'P' => 


"°a 2 tt 


I 


sys 


natural frequency of ballast when attached to moving cone, cps 
natural frequency of cone with nonrigidly mounted ballast, cps 
altitude, ft 

sum of moments of inertia of ballast and cone about transverse 
axes through their centers of gravity, slug-ft 2 

moment of inertia about a transverse axis through the center of 
gravity of the cone-ballast system, slug-ft 2 , 

I + Xv 2 

mi + m2 d 


iii 



k 

Mq 


mi 

m2 


q 

q 

00 

t 

V 

x 

x b 


spring constant, lb/ft 

total pitching moment of cone alone 

mass of cone alone, slugs 

mass of ballast, slugs 

pitching velocity, rad/sec 

free-stream dynamic pressure, lb/ft 2 

flight time, sec 

velocity along flight path, ft/sec 

horizontal displacement from inertial reference, ft 

distance from cone center of gravity (x cg ) to ballast, ft 

& cone 


x c center-of-gravity location of cone alone, ft 

^cone 


Xcg sys 

x cp 

y 

z 

a 

y 

n R 

e 

p 

032 

n 


center-of-gravity location of cone-ballast system, ft 
center-of-pressure location on cone, ft 

displacement of ballast normal to cone axis of symmetry, ft 
vertical displacement from inertial reference, ft 
angle of attack, rad 
swerve angle, rad 

ratio of damping to critical damping for ballast system, 
b 

2/km 2 

pitch angle, rad 
atmospheric density, slugs/ft 3 

natural frequency of cone in equation (13), 2nf , rad/sec 
natural frequency of ballast in equation (13), 2jrf rad/sec 
first derivative with respect to time, 


* 


h 


iv 




(") 


second derivative with respect to time. 


d 2 

dt 2 


( ) £ c frequency coalescence 


4 


I 


V 



THE EFFECTS OF A NONRIGIDLY SUPPORTED BALLAST ON THE 


DYNAMICS OF A SLENDER BODY DESCENDING 
THROUGH THE ATMOSPHERE 
By Gerald N. Malcolm 
Ames Research Center 


SUMMARY 


The effects of a nonrigidly supported ballast on the dynamic behavior of 
a slender body with linear aerodynamics descending through the atmosphere in 
planar motion are investigated. The equations of motion of both the main 
vehicle and the spring-mounted ballast are derived and specific examples are 
presented to illustrate the resulting motions. Divergent oscillations of the 
main vehicle occurred for those conditions near the coalescence of frequency 
of the main body and ballast motions. Divergence occurred earlier in the 
flight when spring stiffness was decreased, or when ballast damping was 
increased for a given spring stiffness. 


INTRODUCTION 


Many atmosphere- entry bodies, particularly slender bodies, are ballasted 
for adequate aerodynamic static stability. The ballast typically is heavy, 
dense, and difficult to support in a perfectly rigid fashion under dynamic con- 
ditions. The purpose of this paper is to report some effects on the vehicle 
dynamics when the ballast is not rigidly supported but oscillates as part of a 
coupled system with the main vehicle. 

The differential equations of planar motion for a body with linear 
aerodynamics and a flexibly mounted mass are derived; and a specific example - 
a slender cone with spring-mass system to represent the ballast arrangement - 
is used to illustrate some of the possible results of vehicle-ballast dynamic 
interaction. 


EQUATIONS OF MOTION 


The system to be analyzed consists of a main body - in this case, a 
slender cone - and an internal nonrigidly supported ballast simulated by a 
simple spring-mass system. The cone motion is restricted to one degree of 
freedom in rotation and two degrees of freedom in translation, all in a single 
plane. The ballast motion is restricted to this same plane and has one- 
dimensional motion perpendicular to the cone axis. The system is shown in 



O Center of pressure 

O Cone center of gravity 

® System center of gravity (with rigid ballast) 




Sketch (b) 


sketch (a) where mi and m 2 are the 
masses of the cone and ballast, 
respectively. The centers of gravity 


of the cone alone, x 
the system, x C g 


and of 


=>sys 


cg cone' 

, (if the ballast 


were rigidly mounted) are shown with 
the aerodynamic center of pressure, 
x C p. The distance of the ballast 


center of gravity from the reference 
point (i.e., the center-of-gravity 
position of the cone alone) is x^. 

Note that x rcr , the center of 
cg sys 

gravity of a system consisting of a 
vehicle with moving parts, does not 
remain stationary so is not a suit- 
able reference. The total spring 
constant in the ballast system is k. 


This system is shown again in 
sketch (b) with an inertial reference 
system and body-fixed aerodynamic 
forces F x and F^. The center of 
rotation for the aerodynamic moment, 
M 0 , is the center of gravity of the 
cone alone. Shown also are angles 
6, the angular position of the cone 
with respect to the inertial x-axis, 
a, the angle of attack, and y, the 
swerve angle. 


Referring to sketch (b) we can write the differential equations of motion. 
If small angles are assumed, for example, 

sin (a) = a 


cos (a) = 1 


r 




and gravity forces are neglected, the moment and force relationships are 


10 = M 0 

+ (ky + by)x b + m 2 xy 

CD 

m l z l = f n 

+ ky + by 

(2) 

1112 Z 2 = ky 

- by 

(3) 

z 2 = Z 1 

+ x b e + y 

(4) 


r 


* 


where 
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NL 


- Cm^q/d + C m q Ad (y) 6 + C md q o!) Ad (v) “ 


F N = C N a “q oo A 


F x = c x q A 


x = - 


mi + m2 


Note that the aerodynamic normal force, F N , and the restoring moment, M 0 , were 
assumed to vary linearly with angle of attack and that the axial force, F x , 
was assumed constant. From sketch (b) 


CD 

II 

P 

+ 

(5) 

0 = a + y 

(6) 

0 = a + y 

(7) 

x = V 

(8) 

• Cxq<»A 

x = V = - 

mi + m 2 

(9) 

zi = Vy 

(10) 

Zi = Vy + Vy 

(11) 

q = (l/2)pV 2 , we solve for y 

and substitute 


for zi from equation (2) 


. _ zi - VY _ C N a apVA ky + b y C x pVA 
Y V 2mi + Vmj + 2 (mi + m 2 ) Y 


( 12 ) 


Differentiating equation (12), we get 


C Nc apVA C Na apVA . C^apVA + b » 

2nii + 2mi + 2m^ Vm! 


+ b y v 


V" 


V + 


C v pVA 


C v pVA 


C„pVA 


2 (mi + m 2 ) 2 (mi + m 2 ) 2 (mi + m 2 ) 


(13) 


If we substitute for V from equation (9) and y from equation (12) and 
neglect second-order terms with products of coefficients such as Cn^C x , 
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.. _ C N a pVA . C N a apVA C x pVA _ k y + by 
Y 2m! a + 2mi + 2(111! + m 2 ) Y + Vmi 


If the terms that are small relative to the other terms in equations (12) and 
(14) are eliminated 

C Na apVA C Na aq oo A 


C Na pVAd Cj^qJVo 


If we now substitute for y and y into equations (6) and (7) and then 
substitute for 0 and 0 in equations (1) and (4), 


c N„q AI 


*■ t l„ Ad (?) (s * c «a) ° 


- C m aq Ad - (ky + by)x, + — - ■ F y = 0 

m ot n oo v J J J b m 1 + m 2 x 


z 2 = zi + x. 


c n “q A \ 

+ -v^-) + y 


Note that the (C m + C m .)a term in equation (15) is damping due to pitching 
and the Cjj a terms in both equations (15) and (16) are damping terms due to 

plunging or swerving of the vehicle. Now equations (2), (3), (15), and (16) 
can be used to produce two simultaneous differential equations in a and y. 
Substituting into equation (16) for and z 2 from equations (2) and (3), 

we get 


? *e=P)&) > * (=^)&) * 


c N„«q m A 


A IN ‘oo A 

+ Xv a + — + x K — r . = 0 

b mi b Vmx 


If we substitute for (ky + by) from equation (17) into equation (15), 


mi + m 2 


A T q A / m n r 

2 ) a - - Cm (I + 1 

I 1N a Vmi \ mi + 


x u 

m 2 b 


+ q-Ad y + c m . j 


- C m aq Ad + — m ^ n ^_ x h y + x k^N + ~ — F x y = ® (18) 

m a » m^ + m 2 b' mi + m 2 b n mi + m 2 x 


Rewriting equations (17) and (18) and substituting coefficients for and F x 
and abbreviating the a coefficient in equation (18) with CMQ, we get 
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* * bp)© * * e=p)&) * 


+ x, a + 


q A 

. ‘no 


-I ■ 

b mi 


C N«(“ * H ■ 0 


( 19 ) 


( : * STTiJ x b 2 )“ - k Ad §)(H“ " ^“(s, • SI ” 2 m 2 - C Na t) a 


m 1^2 


m c 


+ x K y + C v q Ay = 0 (20) 

mi + m 2 o J mi + m 2 xn co J 


Note that in equation (20) {I + [mim 2 /(mi + m 2 )]xb 2 } is the moment of inertia 
about the instantaneous center of gravity of the system, designated I S y S , 
and that {Cm - [m 2 /(mi + m 2 ) ] (x^/d) } is the pitching-moment curve slope 

for the system if the ballast is rigidly mounted, that is, c m as y S - 

In equation (19) b can be expressed conveniently in terms of some 
fraction of critical damping. Critical damping will be arbitrarily defined 
here as that minimum value of b that would cause the ballast to lose its 
oscillatory behavior if it were mounted on a stationary platform. (This is 
not to be construed as the critical damping for the entire system; that quan- 
tity is undefined.) From simple harmonic-motion theory, then, b cr i t = 2/km2 
and 


b = — (2/km 2 ) = nn (2 Am 2 ) 

D crit 


If 2n R A"m 2 is substituted for b in equation (19), then 


m„ + m„ 


y + 2 


mi 


n R J^~2 Y + 


mi 




y + x b“ + ir %(“ + V X b) = ° 


( 21 ) 


Equations (20) and (21) must be solved numerically. However, some 
simplifying assumptions will enable the equations to be solved in closed form, 
and the solutions can be useful in indicating trends. 


Closed-Form Solution 

If q is assumed constant and the aerodynamic and spring damping 
(CMQ = n r “ 0) are neglected, equations (20) and (21) can be rewritten in the 
form 


a + K x a + K 2 y + K 3 y = 0 (22) 

y + K 4 a + K 5 a + K & y = 0 (23) 
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where 


Ki 


-q Ad(C m — — C N --] 

n oo y m a m 1 + m 2 d f 


I + 


mim 


l m 2 


mi + m2 b 


m 


l m 2 


K 2 = 


m i + m 2 ^ 


I + 


n*l m 2 ^ 

Xt 


mi + m2 b 


m,- 


Ks = 


in 1 + m 2 


C x q A 


I + 


m 1 m 2 2 
x b 


mi + m2 


K 4 = 


Cm q A 
"a 1 ® 

mi 


k 5 = *b 


K 6 = 


m l + m 2 k_ 
m l m 2 


Closed-form solutions can be written for these equations, and with initial 
conditions at t = 0, a = a Q , y = y Q , a = a Q , y = y Q they can be expressed as 


y = Ci cos wit + C 2 sin a^t + C3 cos w 2 t + C4 sin w 2 t' 
a = C 5 cos wit + C5 sin wit + C7 cos w 2 t + Cg sin w 2 t 


(24) 


(these expressions are valid only up to and including the point where 
frequency coalescence occurs, that is, a>i = w 2 ) where 


Cl = 


G 


a 


o 


H 


6 


and 


Co = 


r ^ o a o 

b 03 x " Wj 

G - H 


C 3 = y G - Ci 


c 4 = 


y Q - c,oi 


2 1 


0)2 


Cq = ClH 


c 6 = c 2 h 
c 7 = c 3 g 

Cq - C 4 G 


(Kg + - Vh * K s) 

CKi + K 4 ) - (K 5 + l)a ) 2 2 


(Kg + l)t»3i 2 - CK 3 + Kg) 
(K x + K 4 ) - (K s + l )^ 2 


From equation ( 24 ) we can express the primary natural frequencies of the cone 
and ballast, respectively. 


0 )! , 0) 2 


Ai + (A x 2 - A 2 ) 1/2 
2 


(25) 


where 


Ai 


k 3 k 5 + k 2 k 4 - k 6 - k x 

K 2 K 5 - 1 


A 2 = 4 


/k 3 k 4 - k x k 6 \ 

V K 2K 5 - 1 ) 


and K is defined as before. Frequency coalescence is defined as the 
condition where 03 \ = 012- 
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When frequency coalescence takes place, (q /k)^ c can be written in 

terms of the body and aerodynamic- force parameters. From equation (25) it is 
apparent that frequency coalescence occurs (toi = m 2 ) when Ai 2 - A 2 = 0. If 
we substitute the above expressions for Ai and A 2 into A x 2 - A 2 = 0 and 
solve explicitly for (q /k)f c 


Bl 



+ 



= 0 


(26) 


and 


GO* 


-B 2 ± (B 2 2 - 4BiB 3 ) 1/2 


2Bi 


(27) 


where 


Bi = 


m 2 2 A 2 x b 2 C x 2 + 2 A 2 d 2 + " J “ 2 

(mi + m 2) 2 m ct 


2m 2 A 2 dX| ;) 2 C x C ina 4m2A 2 C x Cf,j a I 


mi + m 2 (mi + m2)mi 


Bo = 2A 


B 3 


[c - — m2 ^ + C N ( 2XbI ) + C m d (- 1 + m2 I - x.) 

lx l mi mi + m 2 / 1N a \ mi / m a V mim2 b J 

, f I(mi * * X. 2 1 2 

mim 2 b 


DISCUSSION OF RESULTS 


A complete study of all the parameters that affect the oscillatory motion 
of the vehicle was not attempted. Rather, a specific cone with a given flight 
trajectory, in terms of time, velocity, and dynamic-pressure variation with 
altitude, was used as a model. The effects of spring stiffness and spring 
damping in the simulated ballast system and the effects of aerodynamic damping 
were investigated to ascertain their influence on the oscillatory frequency 
and amplitude of both the ballast and the vehicle. The influence of varying 
the ballast location is also examined briefly. The values used for the geo- 
metric parameters and aerodynamic coefficients are listed below (see 
sketch (a) or (b)): 


d = 2.5 ft 


l = 7.0 ft 

x cg = 3.956 ft 
B sys 
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4.780 ft 


x cg 
6 cone 

x cp = 4.928 ft 

x b = 2.36 ft 

mi = 12.36 slugs 


m 2 = 6.64 slugs 

I = 44.44 slug-ft 2 (I sys = 68.5) 
Cm =1.95 
C x = 0.075 

C m = -0.110 (C m = -0.753) 
a a sys 

^mq + Cm^ = ® > -1.5 


The flight trajectory used for the investigation is shown in figure 1, 
where velocity, dynamic pressure, and oscillatory frequency for a rigid system 
(infinite spring constant) are shown for an altitude range from 200,000 to 
20,000 feet. Also shown is elapsed time with t = 0.0 at h = 200,000 feet. 
The initial entry angle was taken to be -22°. 



Figure 1.- Flight trajectory parameters for example case. 



Oscillatory motions of both conical vehicle and ballast were obtained by 
integrating numerically the coupled equations (20) and (21), where V and q^ 
were entered in tabular form as functions of time. The motions were started 
at the velocity and dynamic pressure corresponding to 200,000 feet altitude. 
The other initial conditions were a = 0.2 rad, a = 0, y = 0, y = 0. The 
resulting motions were then studied to determine the frequency and amplitude 
behavior of both cone and ballast and to determine whether conditions might 
occur at which the cone angle-of-attack oscillations would be divergent. 


Figure 2 shows typical motions of cone and ballast as functions of time. 


For this case k = 5,000 lb/ft, n R = 0, and Cj 


li L aj xl — lx , l/l / v/ O.L//XU, iij^ — L/, cl 1 1 lx 1 '-'in * — vx * Initially the 

primary oscillatory frequencies of the two motions are much different. At 
high altitude (where q^ is very small), the ballast oscillates at a primary 
frequency, f , near the natural frequency for a two-body system uninfluenced 
by any external forces or moment; that is. 


where f 0 is the natural frequency of a simple mass spring-mounted to a 

rigid platform. The frequency is about 6.7 cps, nearly equal to the 6.73 cps 
calculated from equation (28) estimated for flight in a vacuum (where q = 0). 
(The ballast frequency, fcy> would be 4.37 cps if the cone were held station- 
ary.) Initially each motion for y and a has a secondary natural frequency, 
which is simply the other's primary natural frequency. As the dynamic pres- 
sure increases with time, the primary oscillatory frequency, f a , of the cone 
increases at a faster rate than would be expected for a rigid system. The 
difference is shown in figure 3, where f a is plotted as measured directly 
from the motion in figure 2 and f 0 is the frequency for a rigid system. 

Since only a portion of the ballast motion in figure 2 exhibits a primary fre- 
quency that can be measured directly, the fy curve was calculated from the 
closed-form solution (eq. (25)) for various values of q^ and then plotted 
versus the corresponding time, t. The closed-form solution agreed perfectly 
with f a and f that were measurable from figure 2. It is evident that as 
q^ increases (and consequently f a ) , the primary ballast frequency decreases 
until at t = 10.2 seconds the ballast and cone frequencies are the same (fre- 
quency coalescence). However, figure 2 shows that the motion diverges at 
t ~ 9.2 seconds, that is, before frequency coalescence takes place. Despite 
the onset of divergence the frequencies behave as predicted by equation (25). 
The amplitude of both cone and ballast increase until the ballast motion is 
bounded by the cone walls. At this point (t = 10.3 sec, see fig. 2) 
describing the motion of the system becomes more complex and is beyond the 
scope of this investigation. 

Some motions obtained from closed-form solutions help to understand why 
the coupled motion of the ballast and cone cause divergence prior to frequency 
coalescence. Figure 4 illustrates ballast and cone motions at constant values 
of q ro corresponding to values just prior to and at frequency coalescence. 
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Cone frequency with nonrigid ballast (measured from fig. 2) 

- ' Cm °,v S Ad 

0 «. 


Cone frequency with rigid ballast, f^ 
Ballast frequency (eq.(25)) 


isys 



Frequency coalescence (4.87 cps) 



Figure 4.- Oscillatory motion of cone and ballast 
at constant dynamic pressure. 





When the two bodies oscillate at nearly identical natural frequencies 
(fig. 4(a)), a M beat frequency" is established in which the ballast motion 
first lags, then leads the cone motion. At frequency coalescence (fig. 4(b)) 
the ballast lags continuously. Applying this quasi-steady interpretation of 
the lag phenomenon to the variable q^ case in figure 2, we conclude that the 
divergent motion occurs somewhat before frequency coalescence is reached 
because of the >180° lag of the ballast that is typical of part of the beat 
frequency motion. Because q^ is continuously increasing, the beat frequency 
does not fully form and the ballast lags continuously until coalescence. 

The effect of the ballast on the cone motion is illustrated in figure 5 
where the amplitude history of the cone with a nonrigid ballast is compared to 
one calculated for a completely rigid system. Note that for the rigid body 
the amplitude decreases despite the fact that C m ^ + C m ^ = 0 because of the 

increasing q^ with time and the damping contribution from plunging. It is 
apparent that the amplitude ratio is affected by the ballast long before fre- 
quency coalescence conditions occur; in fact, the motion diverges before 
coalescence conditions are reached because of the >180° lag discussed above. 
The point at which the ballast reaches the cone wall is also indicated. 



Figure 5.- Envelope of oscillatory cone motion. 


Effects of Changing Spring Force 
and Damping Constants 

To see what effect a change in spring stiffness or the addition of 
damping to the ballast system has on the cone oscillatory history, amplitude 
ratio (where a Q is 0.2 rad at h = 200,000 ft) versus altitude is plotted in 
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figure 6 for three different spring constants k = 5,000, 15,000, and 
50,000 lb/ft. For comparison, the amplitude ratio for a rigid system is also 
shown. The damping-in-pitch parameter, C m ^ + C m ^, and ballast damping, n^, 

are both zero. The effect of decreasing the stiffness (decreasing k) of the 
ballast is to produce oscillatory divergence at a higher altitude (and 
lower q ) . 

~ co'' 



Figure 6.-. Envelope of oscillatory cone motion for three different spring stiffnesses. 


Ballast damping.- Figure 7 shows the effect of introducing damping into 
the ballast system for three different spring constants. The damping values 
go to a practical upper limit of p R = 1.0 (the critical damping value for the 

ballast system if the ballast were attached to a rigid platform) and the 
effect of increased damping is to cause divergence to occur sooner in the 
flight. The reason is simply that if the ballast system includes damping, the 
ballast will lag behind the cone sooner than it would if there were no 
internal damping. In all cases the amplitude ratio is plotted until the 
ballast contacts the cone wall or the amplitude ratio increases to 1.0, 
whichever occurs first. 

Aerodynamic damping .- Figure 8 shows that aerodynamic damping delays the 
onset of divergence, as one would expect. Note also that for the rigid-body 
case there is no divergence during the decreasing dynamic pressure part of the 
trajectory (i.e., altitudes below 37,000 ft), in either case because even when 
Cm + C m . = 0, there is sufficient damping from the plunging motion (C^j a 

Q Ot CX 

term) to prevent divergence in angle of attack. 


14 






Figure 8.- Effect of aerodynamic damping on cone amplitude history. 


Results From Closed-Form Solutions 


We have already seen in the discussion of figure 3 that the closed-form 
solutions are very useful for calculating frequency histories (eq. (25)) by 
using discrete points along the trajectory and for calculating the point in 
the trajectory where a frequency coalescence condition exists (q can be 

°°fc 

calculated directly from eq. (27) and this in turn can be related to a 
specific time, t, from fig. 1). 


In addition the closed-form solution (eq. (27)) can be used to determine 
the effect of the longitudinal location of the ballast on the coalescence 
point for the system. In figure 9 (q^/k)^ is plotted as a function of the 



x b , ft 

Figure 9.- Variation of dynamic pressure to spring stiffness ratio at frequency coalescence with 

bal last location . 
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ballast location x^. The most important point is that if the ballast 
location is less than x^ = 0,465 ft, there is no theoretical frequency coales- 
cence. Above x^ = 0.465 there are two q^A conditions that will produce 
coalescence although the lower one is more significant since it would be the 
first to be encountered in an atmosphere entry. This coalescence point can, 
of course, be coordinated with a specific altitude by using the q^ history 
with altitude (as in fig. 1). It should be emphasized that this criterion 
for minimum x^ to produce frequency coalescence is for the specific body 
parameters used in this paper and would be numerically different if the body 
dimensions or aerodynamic coefficients in equation (27) were changed. 


CONCLUSIONS 


The effects of a nonrigidly supported ballast on the dynamic behavior of 
a slender body with linear aerodynamics descending through the atmosphere in 
planar motion have been investigated. Some important conclusions can be drawn. 

The point at which frequency coalescence (ballast and cone frequencies 
are equal) occurs in a flight with variable dynamic pressure can be predicted 
accurately by means of closed- form solutions if constant values of dynamic 
pressure are assumed and aerodynamic and ballast damping are ignored. However, 
it has been shown that the amplitude of oscillations diverges before frequency 
coalescence is experienced. This amplitude divergence is caused by the con- 
tinuous >180° phase lag of the ballast motion relative to that of the cone 
just prior to frequency coalescence. (Phase lag is exactly 180° at frequency 
coalescence.) In addition, it was found that 

1. Decreasing the stiffness of the ballast support causes divergence 
earlier in the flight (higher altitude and lower dynamic pressure). 

2. Introducing damping into the ballast system (at least up to the point 
of critical damping for the ballast mounted on a stationary platform) produces 
the phase lag sooner in the flight than if there were no internal damping and 
therefore promotes earlier divergence of the oscillations. 

3. Aerodynamic damping delays the onset of divergence. 

4. Varying the ballast longitudinal location has large effects on the 
altitude at which divergence occurs. There is a minimum distance from the 
cone center of gravity to the ballast location below which no frequency 
coalescence can exist. 
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